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ABSTRACT 

A  direct  numerical  simulation  of  the  Navier-Stokes  equations  of 
a  feedback  controlled  cylinder  wake  flow  at  the  Reynolds  number 
of  60  is  presented.  The  simulation  shows  that  vortex  shedding 
behind  a  circular  cylinder  is  suppressed  to  a  significant  degree 
by  applying  a  feedback  forcing  with  an  appropriate  phase  via 
blowing/suction  slots  located  at  ±110*  from  the  leading  edge. 
Depending  on  the  phase  of  the  feedback  forcing  relative  to  the 
vortex  shedding  phase,  the  shedding  of  the  vortices  is  either 
suppressed  (stabilized)  or  further  enhanced  (destabilized). 

1  INTRODUCTION 

Many  flows  of  engineering  interest  (e.g.  separating  flows)  pro¬ 
duce  the  phenomena  of  vortex  shedding  and  the  associated  struc¬ 
tural  response.  Applications  include  marine  structures,  under¬ 
water  acoustics,  civil  and  wind  engineering  to  name  a  few.  The 
ability  to  tailor  the  wake  of  a  bluff  body  could  be  used  to  reduce 
drag,  increase  mixing  or  heat  transfer  and  enhanced  combus¬ 
tion.  Vortex  shedding  from  bluff  bodies  has  been  reviewed  most 
recently  by  Sarpkaya  (1979)  and  Bearman  (1984).  Vortex  shed¬ 
ding  from  a  cylinder  occurs  over  wide  range  of  Reynolds  number 
from  47  <  Re  <  300,000  and  then  again  at  Re  >  10s.  When 
Reynolds  number  exceeds  a  critical  value  of  approximately  47, 
a  symmetry  breaking  supercritical  Hopf  bifurcation  occurs  lead¬ 
ing  to  a  discrete  frequency  vortex  shedding  known  as  Karman 
vortices.  Chomaz  et  al.  (1988)  have  demonstrated  that  this 
symmetry  breaking  Hopf  bifurcation  is  related  to  the  destabi¬ 
lization  of  a  global  mode.  The  existence  of  an  unstable  global 
mode  is  intimately  linked  to  the  fact  that  the  near  wake  is  abso¬ 
lutely  unstable  followed  by  the  convectively  unstable  far  wake. 
The  understanding  of  the  dynamics  of  two-dimensional  wakes 
at  Reynolds  numbers  in  the  vicinity  of  the  critical  value  is  quite 
firmly  established.  The  preceding  description  of  the  vortex  shed¬ 
ding  mechanism  is  from  an  instability  perspective.  From  a  mech 


anistic  point  of  view,  a  flow-induced  feedback  causes  the  shed¬ 
ding  of  the  vortices,  where  the  changing  location  of  boundary 
layer  separation  on  one  side  of  the  cylinder  causes  a  change  in 
pressure  distribution  that  forces  an  opposite  asymmetric  change 
in  boundary  separation.  Any  small  flow  perturbation  is  enough 
to  cause  the  initial  asymmetry  that  starts  the  feedback  mech¬ 
anism.  The  observation  that  the  vortex  shedding  is  caused  by 
the  moving  location  of  boundary  layer  separation  is  supported 
by  the  absence  of  a  distinct  shedding  frequency  in  the  range 
300,000  <  Re  <  3,000,000  where-'laminar  separation  followed 
by  turbulent  reattachment  fixes  the  turbulent  separation  loca¬ 
tion  so  firmly  that  feedback  is  inhibited. 

Under  the  influence  of  vortex  shedding,  the  changing  pres¬ 
sure  distribution  causes  large  oscillatory  lift  forces.  When  a  cir¬ 
cular  cylinder  is  flexible  and  lightly  damped  or  rigid  and  flexibly 
mounted,  vortex  shedding  can  cause  a  coupled  fluid  structure  in¬ 
teraction.  This  flow-induced  resonance,  causes  the  body  and  the 
wake  to  oscillate  at  the  same  frequency.  This  phenomenon,  al¬ 
ternately  called  galloping,  lock-on  or  phase  locking,  results  in 
spanwise  phase  locking  of  the  fluid  and  the  structure  causing 
coherent  vortex  shedding  along  the  span  of  the  cylinder 

The  important  aspect  of  this  phenomenon,  from  a  control 
standpoint,  is  that  the  large  pressure  forces,  both  lift  and  drag, 
are  a  result  of  the  small  viscous  forces  acting  in  the  bound¬ 
ary  layer  of  the  cylinder  This  implies  that  small  control  forces 
applied  to  the  boundary  layer  can  be  used  to  control  the  larger 
pressure  forces  to  achieve  the  effect  of  modifying  the  vortex  gen¬ 
eration  and  unsteady  lift  forces.  This  concept  has  been  demon¬ 
strated  by  Ffowcs  Williams  ic  Zhao  (1989),  where  the  feedback 
from  a  downstream  hot  film  sensor  was  used  to  drive  a  speaker 
embedded  in  the  wind  tunnel  wall  and  suppress  vortex  shedding 
from  a  cylinder  in  crossflow.  The  acoustic  velocity  perturbations 
from  the  speaker  are  many  limes  smaller  than  the  freest  ream 
velocity,  and  yet  control  was  obtained  This  study  shows  that 
vortex  shedding  control  is  a  real  possibility.  Two  recent  studies 
(Tokomaru  k  Dimolakis  ( 1991 );  Filler  et  al.  (1991 ))  have  shown 
that  small  rotation  oscillations  of  a  circular  cylinder  can  cause 


lock-on.  Numerical  simulation  of  Kanuadakis  Jc  Triantafyllow 
(1!)S9)  have  shown  that  small  variations  of  tin-  near  wake  pro 
(luce  lock-on. 

The  emphasis  in  tins  work  is  on  active  feedback  control,  as 
opposed  to  open  limp  votiUols  as  desculied  in  \\  n  et  ,d  (1989), 
Tokumani  &  Diniolakis  ( 1991 ),  and  Tiller  ct  al.  (1991)  In  these 
studies,  rotary  oscillation  of  the  cyliudei  was  used  to  alter  the 
cylinder  wake,  with  no  feedback  from  the  cyhmlei  or  the  flow 
Open  loop  control  can  cause  lock  on  of  voiles  shedding  and, 
therefore,  modify  tlie  voitcx  shedding  phase,  hut  it  is  highly 
unlikely  that  it  can  suppress  the  vortex  shedding  altogether. 

The  most  powerful  asymmetric  actuator  for  control  of  the 
boundary  layer  of  the  flow  over  a  2-D  cylinder  is  the  rotary  oscil¬ 
lation  of  the  cylinder  wall  itself  (sec  Chang  (1976)  for  numerous 
separation  control  schemes).  However,  this  method  of  boundary 
layer  control  has  certain  practical  problems  in  implementation 
as  well  as  lack  of  generality  to  other  geometries  and  therefore 
will  not  be  considered  in  this  paper. 

The  next  most  powerful  -'actuator  appears  to  be  a  pair  of 
suction/blowing  slots  located  opposite  to  each  other  at  some  lo¬ 
cation  from  the  leading  edge  or  the  stagnation  point.  A  small 
amount  of  fluid  sucked  in  one  slot,  thinning  the  boundary  layer, 
and  ejected  out  through  the  other  slot,  will  asymmetrically  change 
the  location  of  the  boundary  layer  separation. 

In  this  paper,  we  will  consider  the  problem  of  suppressing 
or  stabilizing  the  near  wake  of  a  flow  past  a  circular  cylinder 
at  low  Reynolds  number  by  direct  numerical  simulation  of  2-D 
Navier- Stokes  equations. 


2  NUMERICAL  PROCEDURE 


The  control  problem  is  investigated  by  numerically  solving  the 
2-D  Navier-Stokes  equations  in  a  stream  function/vorticity  for¬ 
mulation  as  follows; 


^  +  v-(^w)  =  |-W  (i) 

vv  =  (2) 

where  w  and  »/»  are  the  vorticity  and  the  stream  function,  respec¬ 
tively.  The  above  equations  are  nondimensionalized  using  the 
cylinder  radius,  R,  the  freestream  velocity,  U,  and  the  kinematic 
viscosity  of  the  fluid,  v.  The  Reynolds  number  Re  is  defined  as 


follows; 


(3) 


where  d  is  the  cylinder  diameter.  The  boundary  conditions  on 
the  cylinder  surface  are 


^(r  =  1,9,0  =  0,  (4) 

0(r=  1,9,0  -  *.(»,*).  (5) 

In  the  case  of  zero  blowing/suction,  0„(fl,f)  is  identically  equal 
to  zero.  Here  we  have  two  boundary  conditions  for  the  stream 
function  while  there  is  no  boundary  condition  for  the  vorticity. 
Therefore,  the  second  order  accurate  boundary  condition  for  the 
vorticity  at  the  cylinder  surface  is  derived  using  equation  (2) 
and  the  boundary  conditions  for  the  stream  function.  The  outer 
boundary  conditions  are 
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advaiicnitcnl  is  tUmc  vi.i  l\vu  h.uUouai  slcps  (Kotc&ak  k  I’ater.i 
(19SG))  The  first  step  is  the  nonhne.tr  explicit  step  to  compute 
the  efTc.c.l  of  the  nonlinear  Atlveciivc  leim  as  follows, 


A/ 


,-0 


(8) 


where  fio  —  23/12,  fit  =  —  1G/12,  and  f}7  =  5/12,  and  tv',  w" 
and  0*  denote  the  intermediate  vorticity,  the  vorticity  and  the 
stream  function  at  n-th  time  level,  respectively.  The  second 
step  is  the  implicit  step,  where  the  viscous  effect  is  computed  as 
follows; 


u>“+1  -w'  2  +u* 

At  ~  Re  2 


(9) 


After  the  vorticity  at  (n+l)-th  lime  level,  w**1,  is  computed,  the 
stream  function,  is  obtained  by  solving  tbe  equation  (2) 


for  tbe  stream  function.  The  accuracy  of  the  numerical  scheme 
is  second  order  in  both  space  and  time. 

The  computational  grid  is  uniformly  stretched  in  radial  di¬ 
rection  such  that  the  grid  is  denser  near  the  cylinder  surface  to 
resolve  the  high  gradient  of  the  flow  there.  A  typical  grid  is 
shown  in  figure  1.  The  outer  computational  boundary  is  located 
at  r»  =  59.38.  The  number  of  grid  points  in  circumferential 
and  radial  directions  are  256  and  128,  respectively. 


3  RESULTS  AND  DISCUSSION 

Accuracy  of  numerical  method  was  validated  by  comparing  the 
numerically  obtained  Strouhal  number  (St)  to  the  experimental 
values  at  the  Reynolds  numbers  of  60,  100,  150  and  200.  The 
nondimensional  frequency  St  is  defined  as  fdd/U,  where  /a  is  the 
frequency  of  vortex  shedding  In  experiments  at  these  Reynolds 
numbers,  due  to  the  finite  aspect  ratio  effect  of  the  cylinder,  the 
vortices  are  shed  from  the  cylinder  at  a  small  oblique  angle  to 
the  cylinder.  Williamson  (1989)  found  that  the  Si  -  Re  dis¬ 
continuities  occurring  al  low  Reynolds  numbers  are  due  to  the 
obliqueness  of  the  vortex  shedding.  He  noticed  that  the  parallel 
shedding  and  the  oblique  shedding  are  related  by  the  cosine  of 
the  angle  as  follows.  The  frequency  of  the  parallel  shedding  can 
be  obtained  by  dividing  the  frequency  of  the  oblique  shedding  by 
the  cosine  of  the  angle  of  the  shedding  that  is  the  angle  between 
the  spanwise  oblique  vortex  line  and  the  cylinder  axis.  He  also 
presented  the  universal  St  -  Re  curve  without  discontinuities  by 
correcting  the  effect  of  the  obliqueness  of  the  vortex  shedding. 
Our  results  are  plotted  with  these  corrected  experimental  val 
ues  of  Williamson  ( 1 989)  in  figure  2.  The  agreement  is  excellent 
The  results  will  presented  for  the  case  of  the  Reynolds  number 
of  60. 

The  lift  coefficient  is  evaluated  directly  from  the  known  vor 
luily  field  l>v  the  following  finimil.i 


60 


CL  =  jr  (  '~7~  +  -~)^i c^(0),w.  (in) 

where  Mi,  is  the  radial  velocity  at  the  cyliiidri  stnf.iii:  due  to 
t|ie  radial  liliiiviug/suction  velocity  tlnnugli  ihe  bfim turn 
slots.  The  first  term  is  the  conti ihulioii  to  the  hit  of  I Fr  l.ui"<  ii- 
toil  viscous  force  and  lh<*  seioud  term  is  the  lifi  dm-  to  the  (lies- 
sure  force.  The  prcsstirr  force  is  the  •loriiniant  lei  in  in  e<|nation 
(10).  In  particular,  f.'(_  is  positive  in  the  positive  <j  duct  lion 

file  local  ioli  of  I  lie  aelnalni  oi  the  lilow  :n^/sm:l  ion  slots 
u  a  critical  factor  in  their  effectiveness  for  the  suppression  of 
the  vortex  shedding.  We  have  tried  the  blowing/siictiou  slots  at 
±50°,  ±05°,  and  ±110”  from  the  leading  edge  of  the  cylinder. 
Witli  t.  c  actuators  at  ±50°,  or  ±65”,  suppression  of  the  vortex 
shedding  was  not  obtained.  We  think  the  ineffectiveness  of  the 
actuator  was  due  to  the  fact  that  the  actuator  was  located  too 
far  from  the  separation  point  on  the  cylinder  surface.  The  sep¬ 
aration  points  on  the  cylinder  surface,  at  the  Reynolds  number 
of  60,  wander  around  ±120*  as  the  vortices  are  shed.  When  the 
blowing/suction  slots  were  located  at  ±110*  from  the  stagna¬ 
tion  point,  that  is  slightly  ahead  of  the  separation  points  of  the 
unforced  wake,  they  proved  to  be  strong  actuators. 

The  blowing/suction  velocity  through  small  slots  is  dictated 
according  to  the  following  feedback  law. 


U>ui 

where  a  is  the  feedback  coefficient,  w(x„y  —  0)  is  the  mea¬ 
sured  vorticity  at  the  downstream  location,  z„  along  the  cen¬ 
terline  and  wm„  is  the  maximum  vorticity  at  the  same  spa¬ 
tial  location  for  the  unforced  wake.  Here,  positive  u/  denotes 
blowing/sucking  through  upper/lower  slots.  This  choice  of  the 
feedback  law  is  motivated  by  the  fact  that  as  the  flow  becomes 
controlled  and  approaches  the  symmetric  configuration,  the  vor¬ 
ticity  u>{z,,y  —  0)  tends  to  zero  resulting  in  the  small  amount 
of  feedback.  Therefore,  as  the  desired  state  is  approached,  the 
required  forcing  level  becomes  small. 

The  phase  relation  of  the  feedback  relative  to  the  vortex 
shedding  phase  was  determined  by  the  downstream  location  of 
the  feedback  “sensor",  x,.  To  properly  choose  x,,  the  move¬ 
ment  of  the  separation  point  as  the  vortices  are  shed  must  be 
known.  Figure  3  is  the  temporal  plot  of  the  vorticity  at  6.01 
radius  downstream  (x  =  -6.01)  along  the  centerline  (y  =  0) 
and  the  vorticity  at  cylinder  surface  (r  =  1)  at  an  angle  of  1 10*, 
which  is  slightly  ahead  of  the  upper  separation  point  for  the  un¬ 
forced  wake.  The  vorticity  measured  on  the  centerline  oscillates 
around  zero  mean  while  the  vorticity  near  the  separation  point 
oscillates  around  some  mean  value.  Let  us  denote  u>c  and  u>,  to 
be  the  vorticity  at  the  centerline  and  the  vorticity  on  the  cylin¬ 
der  surface,  respectively.  The  phase  relation  between  these  two 
are  as  follows; 

ijjc  =  A  sin(2*/(),  (12) 

w,  -  #sin(2*/(  +  <t>)  +  i*w.n,  (13) 

where  d  is  computed  to  be  —0.71 17*.  The  nondimensional  fre¬ 
quency,  /,  is  the  half  of  the  Strouhal  number,  St.  Feeding  back 
the  phase  of  the  vorticity  at  110”  was  tried  first.  The  appropriate 
sensor  location,  z,,  for  this  phase  is  approximately  8  to  9  radius 
downstream.  The  feedback  coefficient  a  was  varied  from  0.3  to 
1.8  for  the  fixed  sensor  location  at  8  5  radius  downstream.  Sup¬ 


pression  of  the  voi  t c x  shedding  was  obtained  although  it  Was  n,,i 
total  suppression  Figure  4  shows  the  vorticity  plots  for  several 
o's.  The  optimum  value  for  a  appears  to  be  close  to  1  I  As  one 
t  leases  from  a  small  value,  the  voiles  shedding  becomes  slightly 
"nppi  e>  -ed,  Milh  maximum  suppression  .Kiitiiinr'  at  I  i  As  o 
becomes  huger  than  II,  lire  voltes  shedding  e,  enhanced.  I  an 
o  I  i.lheneni  w  ,-kr  slimline  looks  mo  ,1  -.  \  muu  t  r  k  al  Fig 
•  Mr’  *1  .Or-  lire  i  an  lesp, nulling  si  le.iinlme  pl,,l>  I'm  the  nplimum 
rail:,  of  u.  I  hr  ter  in  Iil.il  mg  /one  m  elongated  appio.u  hmc  llie 
sle.idr  solution  A  i  •  oi  ding  to  ( In-  i  ale  idat  Ion  .if  lol  itbei  g  (  |  'JtfOj 
of  steady  llmv  past  a  circulai  cylinder  at  He  Ml,  the  length 
of  the  recirculating  wake  bubble  is  around  4  cylinder  diameters 
When  the  shedding  is  nearly  suppressed,  the  recirculating  wake 
(nibble  is  approximately  3  to  4  cylinder  diameter  long 

Next,  we  fixed  the  feedback  coefficient  o  at  I  !  and  varied 
the  sensor  location  The  suppression  of  vortex  shedding  was 
obtained  only  at  narrow  range  of  the  sensor  location,  which  is 
from  7.5  radius  downstream  to  9  radius  downstream.  The  opti¬ 
mum  sensor  location  is  around  8  radius  downstream.  In  figure  6, 
the  vorticity  plots  for  the  flows  with  the  sensors  located  at  vari¬ 
ous  locations  are  given.  The  corresponding  streamline  plots  are 
shown  in  figure  7.  With  the  sensor  located  at  two  radius  down¬ 
stream,  the  vortex  shedding  process  is  enhanced  instead  of  being 
suppressed  as  shown  in  figure  6a.  As  the  forcing  is  turned  on, 
the  vorticity  and  velocity  fluctuations  increase  by  about  100% 
within  one  vortex  shedding  period.  Ffowcs  Williams  ti  Zhao 
(1989)  also  observed  that  the  phase  reversal  of  the  loud  speaker 
from  the  controlled  case  brought  about  an  enhancement  of  the 
fluctuation.  Let  us  take  a  look  at  the  optimum  case,  which  is 
the  flow  with  the  feedback  sensor  located  at  4  diameters  down¬ 
stream.  The  near  wake  vorticity  plot  in  figure  6b  looks  almost 
symmetric.  The  streamline  plot  in  figure  7b  shows  the  recir¬ 
culating  wake  bubble  of  the  length  of  approximately  4  cylinder 
diameters.  1-igure  8  and  9  are  the  time  series  of  the  lift  coeffi¬ 
cient  and  the  streamwise  velocity  at  one  diameter  downstream 
and  one  radius  above  the  centerline,  respectively  Forcing  is 

turned  on  at  t  =  34  and  is  turned  off  at  t  =  500.  As  the  forc¬ 
ing  is  turned  on,  the  fluctuations  decay  rapidly  to  a  low  level. 
The  lift  coefficient  fluctuation  amplitude  and  the  velocity  fluc¬ 
tuation  amplitude  are  reduced  to  25%  and  15%  of  the  unforced 
case,  respectively.  The  mean  of  the  fluctuating  velocity  has  also 
decreased.  As  forcing  is  turned  on,  the  cylin  der  wake  slowly  re¬ 
turns  to  the  natural  vortex  shedding  state  as  also  observed  by 
Ffowcs  Williams  it  Zhao  (1989).  Figure  10  is  the  time  series  plot 
of  the  square  of  the  lift  coefficient  on  a  semi-log  scale  As  forcing 
is  turned  off,  the  fluctuation  begins  to  increases  exponentially  as 
demonstrated  by  the  straight  line  on  a  semi-log  plot  This  region 
is  called  a  linear  growth  region  which  is  followed  by  a  transition 
region  and  a  saturation  state  as  also  observed  in  the  numerical 
simulation  by  Hannemann  Se  Oertel  (1989).  When  the  forcing  is 
on,  the  fluctuations  modulate.  The  modulations  of  the  fluctua 
tions  can  be  explained  as  follows.  As  the  forcing  becomes  small, 
the  suppression  of  the  vortex  shedding  relaxes  slightly  resulting 
in  the  increased  vorticity  fluctuation.  This  increases  the  forcing 
level  once  again  The  increased  forcing  suppresses  the  vortex 
shedding  process.  Then  the  forcing  decreases  again 

Lastly,  we  briefly  look  at  the  downstream  development  of  the 
vortical  structures  Figure  1 1  shows  the  vortical  development  of 
the  wakes  of  the  unforced  case  and  that  of  the  suppressed  case 
up  to  27  diameters  downstream  Although  the  near  wake  is  sup 
pressed  to  a  significant  degree,  the  small  fluctuation  in  the  neai 
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w  wake  becomes  amplified  due  to  the  convective  nature  of  the  in¬ 
stability  of  the  far  wake  as  the  fluctuation  travels  downstream. 
This  is  also  observed  by  Ffowes  Williams  Zhao  (19S9)  and 
shown  in  figure  5  of  these  authors.  Despite  of  the  fact  that  the 
near  wake  spectral  level  is  low  si  lire  the  wake  is  controlled,  the 
spectral  level  increases  exponentially  in  the  downstream  direc¬ 
tion. 

4  CONCLUSION 

We  have  been  able  to  suppress  the  vortex  shedding  behind  a 
circular  cylinder  in  a  two  dimensional  numerical  simulation  at 
the  Reynolds  number  of  60  by  active  feedback  control  method. 
The  actuators  we  chose  to  use  were  the  blowing/suction  slots 
located  at  ±110“*  from  the  leading  edge  of  the  cylinder.  The 
control  was  possible  in  only  a  narrow  range  of  phase  while  the 
vortex  shedding  was  enhanced  if  forced  with  the  phase  out  of 
this  range.  Although  the  present  results  are  for  the  Reynolds 
number  of  60,  a  quite  low  Reynolds  number,  the  authors  be¬ 
lieve  that  this  active  feedback  control  will  still  be  operative  even 
at  higher  Reynolds  numbers.  The  higher  Reynolds  number  flow 
control  is  currently  being  studied.  One  of  us  (D.  S.  Park)  thanks 
the  support  of  the  ONT  (Office  of  Naval  Technology)  Postdoc¬ 
toral  Fellowship. 
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Figure  1:  Computational  grid;  256  x  128. 


Figure  2:  St  versus  Re,  solid  line  is  data  from  Williamson  (1989) 
and  solid  square  is  the  present  data 
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Figure  4:  Vorticity  contour  for  the  sensor  (x,)  located  at  8.5 
radius  downstream  (*,  =  —8.5).  Flow  is  from  right  to  left,  a) 
o  =  0  (natural  shedding),  b)  a  =  0.3,  c)  a  =  1.1,  d)  a  =  1.8. 


Figure  6:  Vorticity  contour  for  o 
c)  x,  =  -9. 


(c) 


Figure  7:  Streamline  plot  for  a  =  1.1.  a)  x,  —  —2,  b)  x,  =  -8, 
c)  x,  =  -9. 


Figure  10:  Scui-log  plot  <rf  time  series  of  lire  square  of  lift  co¬ 
efficient  for  a  =  1 . 1  and  x,  =  8.  The  response  after  forcing  is 
turned  off  is  shown. 


(a) 


Figure  8:  Plot  of  time  series  of  lift  coefficient  for  a  =  1.1  and 
x,  —  —8.  Forcing  is  turned  on  at  t  =  34  and  turned  off  at 
t  =s  500. 
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Figure  9:  Plot  of  time  series  of  Streamwisc  velocity  at  x  =  -2 
and  y  =  I  for  o  =  1.1  and  x,  —  —8.  Forcing  is  turned  on  at 
t  -  34  and  turned  off  at  t  =  500. 


Figure  11:  Downstream  development  of  vorticily.  a)  natural 
shedding,  b)  suppressed  shedding  (or  =  1.1,  x,  =  —8). 


